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We calculate the probability distribution of the transmission eigenvalues T n of Bogoliubov quasi- 
particles at the Fermi level in an ensemble of chaotic Andreev quantum dots. The four Altland- 
Zirnbauer symmetry classes (determined by the presence or absence of time-reversal and spin- 
rotation symmetry) give rise to four circular ensembles of scattering matrices. We determine P({T n }) 
for each ensemble, characterized by two symmetry indices /3 and 7. For a single d-fold degenerate 
transmission channel we thus obtain the distribution P{g) oc g -1+/3/ ^(l — g) 1 ^ 2 of the thermal con- 
ductance g (in units of dir 2 k%To/6h at low temperatures To). We show how this single-channel limit 
can be reached using a topological insulator or superconductor, without running into the problem 
of fermion doubling. 

PACS numbers: 74.25.fc, 05.4S.Mt, 65.80.-g, 74.45. +c 



I. INTRODUCTION 

The Landauer approach to quantum transport^— re- 
lates a transport property (such as the electrical or ther- 
mal conductance) to the eigenvalues T n of the transmis- 
sion matrix product tt'. If transport takes place through 
a region with chaotic scattering (typically a quantum 
dot), random- matrix theory (RMT) provides a statis- 
tical description. 4-6 While the properties of individual 
chaotic systems are highly sensitive to the microscopic 
parameters of the scattering region, such as its geometry 
or the arrangements of impurities, they obey universal 
statistical features, independent of these details, on en- 
ergy scales below the Thouless energy (the inverse of the 
dwell time). The distribution P({T n }) of the transmis- 
sion eigenvalues then naturally emerges as the determin- 
ing quantity for the distribution of the transport proper- 
ties. 

While microscopic details do not influence the statis- 
tics, the role of symmetries is essential. According to 
Dyson^^ there are three symmetry classes in normal 
(non-superconducting) electronic systems, characterized 
by a symmetry index /? depending on the presence or 
absence of time-reversal and spin-rotation symmetry (cf. 
Table [XJ) . The transmission eigenvalue distribution for 
these three RMT ensembles is known^iiS For a single d- 
fold degenerate channel at the entrance and exit of the 
quantum dot this gives the distribution 

P(g) cx g~ 1+ ^ 2 , 0<g<l, (1) 

of the electrical conductance g (in units of de 2 /h). The 
full distribution P({T n }) has found a variety of physical 
applications jii and has also been used in a more mathe- 
matical context to obtain exact results for electrical con- 
ductance and shot nois o 12 : 13 and to uncover connections 
between quantum chaos and integrable models^ 

As first shown by Altland and Zirnbauer^ Dyson's 
classification scheme becomes insufficient in the presence 
of superconducting order: The particle-hole symmetry 
of the Bogoliubov-De Gennes Hamiltonian produces four 



new symmetry classes*^ - — Depending again on the pres- 
ence or absence of time-reversal and spin-rotation sym- 
metry, these classes are characterized by /3 and a second 
symmetry index 7 (cf. Table HI)) . 19,2 ° As we show in this 
paper, the analogous result to Eq. ([T]) is 

P(g)K 9 - 1+ ^ 2 (l~gT/ 2 , 0<.g<l, (2) 

where now g is the thermal conductance in units of 
dn 2 kgT /6h (at temperature To)- We consider thermal 
transport instead of electrical transport because the Bo- 
goliubov quasiparticles that are transmitted through a 
superconducting quantum dot carry a definite amount of 
energy rather than a definite amount of charge. (Charge 
is not conserved upon Andreev reflection at the super- 
conductor, when charge-2e Cooper pairs are absorbed by 
the superconducting condensate.) 

Concerning previous related studies, we note that the 
electrical conductance has been investigated by Altland 
and Zirnbauer j 1 ^ but not the thermal conductance. Ther- 
mal transport in superconductors has been studied in 
connection with the thermal quantum Hall effect in 
two dimensions and also in connection with onc- 
dimensional localization^ 4 * 2 ^ The present study comple- 
ments these works by addressing the zero-dimensional 
regime in connection with chaotic scattering. 

The outline of this paper is as follows. Sections [TT1 
and IIIII formulate the problem and present P({T n }). In 
Sec. IIVI we then apply this to the statistics of the ther- 
mal conductance. The probability distribution ([2]) in the 
single-channel limit is of particular interest (since it is 
furthest from a Gaussian), but it can only be reached in 
the Andreev quantum dot in the presence of spin-rotation 
symmetry. A fermion-doubling problem stands as an ob- 
stacle when spin-rotation symmetry is broken. We show 
how to overcome this obstacle in Sec. |V] using topolog- 
ical phases of matter— ~— (topological superconductors 
or insulators). We close in Sec. I VII with a summary and 
a proposal to realize the superconducting ensembles in 
graphene. 
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Ensemble name 


CUE 


COE 


CSE 


Symmetry class 


A 


AI 


All 


S-matrix elements 


complex 


complex 


complex 


S-matrix space 


unitary 


unitary symmetric 


unitary selfdual 


Time-reversal symmetry 


X 


S = S T 


S = CT25 T (J2 


Spin-rotation symmetry 


x or / 


J 


X 


degeneracy d of T n 


1 or 2 


2 


2 




2 


1 


4 



TABLE I: Classification of the Wigner-Dyson scattering matrix ensembles for normal (non-superconducting) systems, with the 
parameter ft in the distribution flj of the electrical conductance. (The parameter 7 = in these ensembles.) The abbreviations 
C(U,0,S)E signify Circular (Unitary,Orthogonal,Symplectic) Ensemble. The Pauli matrix aj acts on the spin degree of freedom. 



Ensemble name 


CRE 


T-CRE 


CQE 


T-CQE 


Symmetry class 


D 


Dili 


C 


CI 


S-matrix elements 


real 


real 


quaternion 


quaternion 


S-matrix space 


orthogonal 


orthogonal selfdual 


symplectic 


symplectic symmetric 


Particle-hole symmetry 


S = S* 


S = S* 


S = T2S*T2 


S — T2S*T2 


Time-reversal symmetry 


X 


S — (T2S T Cr2 


X 


S = S T 


Spin-rotation symmetry 


X 


X 


/ 




degeneracy d of T n 


1 


2 


4 


4 


P 


1 


2 


4 


2 


7 


-1 


-1 


2 


1 



TABLE II: Classification of the Altland-Zirnbauer scattering matrix ensembles for superconducting systems. For each ensemble 
the parameters f3, 7 in the distribution ([2]) of the thermal conductance are indicated. The Pauli matrices Oj and Tj act on, 
respectively the spin and particle-hole degrees of freedom. The abbreviations (T)-C(R,Q)E signify (Time-reversal-symmetric)- 
Circular (Real, Quaternion) Ensemble. 

tors. (See Ref. 30 for a review.) We assume s-wave su- 
perconductors, with an isotropic gap A, so for excitation 
energies E < A there are no modes propagating into the 
superconductors. In order to enable quasiparticle trans- 
port, the cavity has two additional leads connected to it 
which support Ni , N2 propagating modes (not counting 
degeneracies). The leads connect the cavity to normal- 
metal reservoirs in local thermal equilibrium. 

Quasiparticle transmission is possible only if the exci- 
tations of the Andreev quantum dot (without the leads) 
are gapless. This is also necessary for the excitations to 
explore the phase space of the cavity, an essential require- 
ment for chaotic scattering. Gapless excitations are en- 
sured by taking two superconducting electrodes with the 
same contact resistance and a phase difference tt. This 
value of the phase difference closes the gap while respect- 
ing time-reversal invariance (because phase differences tt 
and —tt are equivalent). Time- reversal invariance can be 
broken by application of a magnetic field, perpendicular 
to the plane of the dot. (A sufficiently strong magnetic 
field closes the gap, so then the 7r-phase difference of the 
superconductors is not needed and a single superconduct- 
ing electrode is sufficient.) Spin-rotation symmetry can 
be broken by spin-orbit coupling. An ensemble of chaotic 
systems can be generated, for example, by varying the 




FIG. 1: Quantum dot in a two-dimensional electron gas, con- 
nected to a pair of superconductors (shaded) and to two 
normal-metal reservoirs. One of the normal reservoirs is at 
a slightly elevated temperature To + ST. 

II. FORMULATION OF THE PROBLEM 

A. Andreev quantum dot 

An Andreev quantum dot, or Andreev billiard, is a 
confined region in a two-dimensional electron gas con- 
nected to superconducting electrodes (see Fig. [1]). Elec- 
tronic transport through this system is governed by the 
interplay of chaotic scattering at the boundaries of the 
quantum dot and Andreev reflection at the superconduc- 
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shape of the quantum dot or by a random arrangement 
of impurities. 

In global equilibrium the superconducting and normal- 
metal contacts are all at the same temperature Tq and 
Fermi energy (or chemical potential) Ep. For thermal 
conduction in the linear response regime we raise the 
temperature of one of the normal metals by an amount 
ST <C Tq . The thermal conductance G is the heat current 
between the normal reservoirs divided by ST. (The reser- 
voirs are kept at the same chemical potential, so there is 
no thermo-electric contribution to the heat current.) 

If fceXb is small compared to the Thouless energy (the 
inverse dwell time in the quantum dot), then G is de- 
termined by the transmission eigenvalues at the Fermi 
energy, 

G = dG J2 T n- (3) 

n 

The sum runs over the min (N\, N 2 ) nonzero transmis- 
sion eigenvalues T n , with spin and/or particle- hole de- 
generacy accounted for by the factor d. The thermal 
conductance quantum for superconducting systems is 
Go = Tr 2 kgTo/6h, one-half the normal-state value^^ 

B. Scattering matrix ensembles 

The scattering matrix 5 is a unitary matrix of dimen- 
sion (Ni + N 2 ) x (Ni + N 2 ) that relates the amplitudes 
of outgoing and incoming modes in the two leads con- 
nected to the normal reservoirs. The energy is fixed at 
the Fermi level (E = 0). Four sub-blocks of 5 define the 
transmission and reflection matrices, 

5=| r ' Nl x Nl ^' Ni xN ?\ (4) 

\tN2XNi r 'N 2 xN 2 J 

(The subscripts refer to the dimension of the blocks.) 
Table [TT| lists the Altland-Zirnbauer symmetry classes 
to which 5 belongs, and the corresponding RMT 
ensembles.— ~— We briefly discuss the various entries in 
that table. 

In the case of systems without spin-rotation sym- 
metry, it is convenient to choose the Majorana basis 
in which 5 has real matrix elements^ Without time- 
reversal symmetry (symmetry class D), the scattering 
matrix space is thus the orthogonal group. The pres- 
ence of time-reversal symmetry imposes the additional 
constraint 5 = (i2S T <J2, where o~j is a Pauli matrix in 
spin-space, and T indicates the matrix transpose. The 
scattering matrices in this symmetry class Dili are self- 
dual orthogonal matrices. (The combination o~ 2 A T a 2 is 
the so-called dual of the matrix A.) 

If spin-rotation symmetry is preserved, the spin degree 
of freedom can be omitted if we use the electron-hole ba- 
sis (rather than the Majorana basis). The electron- hole 
symmetry relation then reads 5 = T25*T2, where now 



the Pauli matrices tj act on the electron-hole degree of 
freedom. The matrix elements of 5 can be written in 
the quaternion form aoTo + i X} n =i a n T m with real coeffi- 
cients a n . The scattering matrix space for the symmetry 
class C without time-reversal symmetry is the symplec- 
tic group, additionally restricted to symmetric matrices 
in the presence of time- reversal symmetry (class CI). 

Henceforth we assume that the quantum dot is con- 
nected to the leads via ballistic point contacts. The RMT 
ensembles in this case are defined by 5 being uniformly 
distributed with respect to the invariant measure dfi(S) 
in the scattering matrix space for each particular sym- 
metry classJ£ (For the distribution in the case that the 
contacts contain tunnel barriers, see Ref. I32I .) 

It is convenient to have names for the Altland- 
Zirnbauer ensembles, analogous to the existing names for 
the Dyson ensembles. Zirnbauer— has stressed that the 
names D, Dili, C, CI given to the symmetry classes (de- 
rived from Cartan's classification of symmetric spaces) 
should be kept distinct from the ensembles, because a 
single symmetry class can produce different ensembles. 
Following Ref. [33, we will refer to the Circular Real En- 
semble (CRE) and Circular Quaternion Ensemble (CQE) 
of uniformly distributed real or quaternion unitary matri- 
ces. The presence of time-reversal symmetry is indicated 
by T-CRE and T-CQE. (The prefix T can also be thought 
of as referring to the matrix transpose in the restrictions 
imposed by time- reversal symmetry.) 

III. TRANSMISSION EIGENVALUE 
DISTRIBUTION 

A. Joint probability distribution 

Because of unitarity, the matrix products ttf and t't 
have the same set T\, T2, . . . T^ miti of nonzero eigenval- 
ues, with iV m i n = min (JVi, -/V2). The calculation of the 
joint probability distribution P({T n }) of these transmis- 
sion eigenvalues from the invariant measure d/i(5) is out- 
lined in App. lAl^ (It is equivalent to the calculation of 
the Jacobian given in Ref. |2J.) The result is 

P({T n }) oc Y[T^ /2)(Nl - N2) Tr 1+ ^ 2 (l - T ( ) 7/2 

i 

xni r *- r /- ( 5 ) 

j<k 

The values of the parameters f3 and 7 characterizing the 
Altland-Zirnbauer symmetry classes are listed in Table 

The distribution ([5]) differs from the result^^ in the 
Dyson ensembles by the factor Y[%(^ ™?i) 7 ^ 2 - Depending 
on the sign of 7, this factor produces a repulsion or at- 
traction of the Ti's to perfect transmission. In contrast, 
the factor T[ i y~ 1+,3/2 j wmc h exists also in the Dyson en- 
sembles, repels or attracts the Tj's to perfect reflection. 
The distributions P{T\) for N\ — N 2 = 1 in the various 
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FIG. 2: Probability distribution © in the case Ni = N 2 = 1 
of a single (d-fold degenerate) transmission eigenvalue T, 
which then corresponds to the (dimensionless) thermal con- 
ductance g — G/dGo- The four curves correspond to the four 
superconducting ensembles in Table [TT1 




FIG. 3: Transmission eigenvalue densities in the T-CQE for 
various numbers N = Ni — N2 of transmission eigenvalues, 
calculated from Eq. ([5]). The large- TV limit is the same for 
each ensemble. 



ensembles are plotted in Fig. [2] In view of Eq. this 
is just the distribution ([2]) of the thermal conductance in 
the single-channel limit announced in the Introduction. 
(How to actually reach this limit is discussed in following 
Sections.) 



B. Eigenvalue density 

The density p(T) of the transmission eigenvalues is de- 
fined by 



(6) 



where (■ • • ) denotes an average with distribution ([5]). It 
can be calculated for Ni , N 2 ^> 1 using the general meth- 
ods of RMT4 

To leading order in Ni , N 2 the eigenvalue density ap- 
proaches the j3 and 7 independent limiting form^&iS 



Po(T) 



Ni + N 2 (T - T ( 



1/2 



2tt \1-T J T 
x 9(1 - T)6(T - T c ), 
{Ni - N 2 ) 2 
(JVi+JVb) 8- 



(7) 
(8) 



(The function Q(x) is the unit step function, Q(x) = 
if x < and Q(x) = 1 if x > 0.) The approach to this 
ensemble-independent density with increasing Ni = N 2 
is shown in Fig. [3] for one of the ensembles. 

The first correction Sp to p is of order unity in Ni,N 2 , 
given by 



Sp(T) 



lO.-2/0)[8(l-T)-8(T-T o )] 
- 1(^/(3)6(1 -T) 

1 e(i-r)e ( r- re ) 

2ir y,!HJ ^(l-T)(T-T c ) 



(9) 



We will use this expression in Sec. HVBI to calculate the 
weak localization effect on the thermal conductance. 



IV. DISTRIBUTION OF THE THERMAL 
CONDUCTANCE 

A. Minimal channel number 

The strikingly different probability distributions ([TJ 
and ([2]) in the normal and superconducting ensembles ap- 
ply to transmission between contacts with a single (pos- 
sibly degenerate) non- vanishing transmission eigenvalue. 
For the normal ensembles a narrow point contact suffices 
to reach this single-channel limit. In the superconducting 
ensembles a narrow point contact is not in general suf- 
ficient, because electrons and holes may still contribute 
independently to the thermal conductance. 

Consider the Andreev quantum dot of Fig. [TJ The 
minimal number of propagating modes incident on the 
quantum dot from each of the two leads is 2 x 2 = 4: 
a factor-of-two counts the spin directions, and another 
factor-of-two the electron-hole degrees of freedom. In 
the CQE and T-CQE the four transmission eigenvalues 
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correction. From Eqs. (JT))-© we obtain 



FIG. 4: Probability distribution of the dimensionless ther- 
mal conductance in the two ensembles with broken spin- 
rotation symmetry, for two independent transmission eigen- 
values (TVi = A?2 = 2). This is the minimal channel number 
in an Andreev quantum dot. To reach the single-channel case 
in the CRE or T-CRE (Ni = N 2 = 1, plotted in Fig. [2} one 
needs a topological phase of matter, as discussed in Sec. [V] 



are all degenerate, so we have reached the single-channel 
limit where the distribution ^) applies. 

The situation is different in the CRE and T-CRE. 
In the T-CRE two of the four transmission eigenvalues 
are independent (and a two-fold Kramers degeneracy re- 
mains). In the CRE all four transmission eigenvalues 
are independent, but two of the four can be eliminated 
by spin-polarizing the leads by means of a sufficiently 
strong magnetic field. So the case with two independent 
transmission eigenvalues (with degeneracy factor d = 2 
for the T-CRE) is minimal in the Andreev quantum dot 
with broken spin-rotation symmetry. 

We have calculated the corresponding probability dis- 
tribution of the (dimensionless) thermal conductance 
g = T\ + T 2 by integrating over the transmission eigen- 
value distribution ([3]). The result, plotted in Fig. [U has 
a singularity at g = 1, in the form of a divergence in the 
CRE and a cusp in the T-CRE. It is entirely different 
from the distribution in the single- channel case (see Fig. 
[5]). How to reach the single-channel limit in the CRE and 
T-CRE using topological phases of matter is described in 
Sec. m 



B. Large number of channels 

In the limit N% , N2 3> 1 of a large number of channels 
the distribution of the thermal conductance is a narrow 
Gaussian. We consider first the average and then the 
variance of this distribution. 

The average conductance can be calculated by inte- 
grating over the eigenvalue density p(T) of Sec. IIIIBI 
We write the average of the dimensionless thermal con- 
ductance g — G/dGo as (g) = go + Sg, where go is the 
leading order term for large N% , N2 and Sg is the first 



9o = 



N X N 2 



N 2 



1 , . 7ViiV 2 

S9= ^- 2 - j) (N 1+N2 r 



(10) 

(11) 



The result (1111) for Sg in the zero-dimensional regime of a 
quantum dot has the same dependence on the symmetry 
indices as in the one-dimensional wire geometry studied 
by Brouwer et al<2i 

Filling in the values of /3, 7, and d in the four su- 
perconducting ensembles from Table lU we see that (for 
Ni = N 2 ) 



in the CRE and CQE, 

SG = <( -Go/2 in the T-CQE, 
G /4 in the T-CRE. 



(12) 



This is fully analogous to the weak (anti)localization ef- 
fect for the electrical conductance (with Go = e 2 /h) 
in the non-superconducting ensembles. 4 Without time- 
reversal symmetry (in the CRE, CQE, and CUE) there 
is no effect (SG — 0), with both time-reversal and spin- 
rotation symmetry (in the T-CQE and COE) there is 
weak localization (SG < 0) and with time-reversal sym- 
metry but no spin-rotation symmetry (in the T-CRE and 
CSE) there is weak antilocalization (SG > 0). 

Turning now to the variance, we address the thermal 
analogue of universal conductance fluctuations. It is a 
central result of RMT 4 that the Gaussian distribution 
of g has a variance of order unity in the large TV- limit, 
determined entirely by the eigenvalue repulsion factor 
lli<j iTi-Tjf in the probability distribution ©. The 7- 
dependent factors plays no role. The result of the Dyson 
ensembles AiS 



Varg 



2(A^ 2 ) 2 
P(Ni + N 2 



(13) 



therefore still applies in the Altland-Zirnbauer ensembles. 

For A^x = N 2 we find the variance of the thermal 
conductance VarG = G^/p with p = 8,4,2,1 in, re- 
spectively, the CRE, T-CRE, CQE, T-CQE. Breaking of 
time-reversal symmetry thus reduces the variance of the 
thermal conductance in the superconducting ensembles 
by a factor of two, while breaking of spin-rotation symme- 
try reduces it by a factor of four. This is fully analogous 
to the electrical conductance in the non-superconducting 
ensembles. 



C. Arbitrary number of channels 

While the results from the previous subsection for the 
average and variance of the thermal conductance hold in 
the limit of a large number of channels, it is also possible 
to derive exact results for arbitrary Ni,N 2 . Following 
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Px + ^Py 



Px ~ iPy 




FIG. 5: Realization of single-channel transmission in the 
CRE, following Ref. l33l . The arrows indicate the direction 
of propagation of chiral Majorana modes at the edges of a 
Px ± ipy-wecve superconductor. The shaded strip at the cen- 
ter represents a disordered boundary between two domains 
of opposite chirality. The thermal conductance is measured 
between two reservoirs at a temperature difference ST, and 
has the single-channel distribution ((2)1 (with ft = 1, 7 = —1). 



FIG. 6: Realization of single-channel transmission in the T- 
CRE. The conducting surface of a topological insulator is 
partially covered by an s-wave superconductor, with order 
parameter ±Ao. Two contacts at temperature difference ST 
inject quasiparticles via two pairs of helical Majorana modes 
(indicated by arrows). For chaotic scattering in the central 
region, the thermal conductance is given by the single-channel 
distribution © (with /3 = 2, 7 = -1). 



the method described in Ref. [TJ, the moments of g can 
be evaluated using the Selberg integral^ We find 



<<?> = 



NiN2_ 



Var g 



2JViJV 2 (iVi+0(iV2+0 
p{N t - 1 + $(N t + t) 2 (N t + £ + 2 IP) ' 



(14) 
(15) 



where we abbreviated N t = N1+N2 and £ = (2-/3+7) / P- 
One readily checks that the large- N limits (|TU| . (jTTJ) , and 
(fj"3f are consistent with Eqs. (M and (fi"5|) . 



HOW TO REACH THE SINGLE-CHANNEL 
LIMIT USING TOPOLOGICAL PHASES 



As explained in Sec. IIV A[ the single-channel distribu- 
tion @ of the thermal conductance can only be realized 
in an Andreev quantum dot in two of the four super- 
conducting ensembles: CQE and T-CQE. The minimal 
channel number in the CRE and T-CRE is two, with 
an entirely different conductance distribution (compare 
Figs. [2] and |4}. Here we show how this fermion doubling 
can be avoided using topological insulators or supercon- 
ductors. 

Consider first the CRE. To have just a single nonzero 
transmission eigenvalue we need incoming and outgo- 
ing modes that contain only half the degrees of freedom 
of spin-polarized electrons. These so-called Majorana 
modes propagate along the edge of a two-dimensional 
spin-polarizcd-triplet, p x ±ip y -w&ve superconductori 26 i 34 
Following Ref. 33|, we consider the scattering geometry 
shown in Fig. [5] The role of the quantum dot is played 
by a disordered domain wall between p-wave supercon- 
ductors of opposite chirality. The system has two in- 
coming and two outgoing Majorana modes, with a 2 x 2 
scattering matrix in the CRE. The thermal conductance 



between the two domains has the single-channel distri- 
bution © (with f} = 1, 7 = -1). 

We now turn to the T-CRE. For a single two-fold de- 
generate transmission eigenvalue we need a 4 x 4 scat- 
tering matrix. Time-reversal invariant scattering in this 
single-channel limit can be achieved if one uses helical 
Majorana modes (propagating in both directions) instead 
of chiral Majorana modes (propagating in a single direc- 
tion only). These can be realized using s-wave super- 
conductors deposited on the two-dimensional conducting 
surface of a three-dimensional topological insulator^ 

The scattering geometry is illustrated in Fig. [51 The 
helical Majorana modes propagate along a channel with 
superconducting boundaries having a phase difference of 
7r (order parameter ±Ao). Two normal-metal contacts 
at a temperature difference ST inject quasiparticles via 
a pair of these modes into a region with chaotic scat- 
tering (provided by irregularly shaped boundaries or by 
disorder) . The ir phase difference of the superconductors 
that form the boundaries of the quantum dot also ensures 
that there is no excitation gap in that region. There are 
four incoming and four outgoing Majorana modes, so the 
scattering matrix has dimension 4x4 and the thermal 
conductance has the single-channel T-CRE distribution 
© (with p = 2, 7 = -1). 

The geometry of Fig.[6]also provides an alternative way 
to reach the single-channel limit in the CRE. One then 
needs to replace the two superconducting islands having 
order parameter — Ao by ferromagnetic insulators. The 
Majorana modes transform from helical to chiral^ and 
one has essentially the same scattering geometry as in 
Fig. [S] — but with s-wave rather than p-wave supercon- 
ductors. 
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VI. CONCLUSION 

In conclusion, we have obtained the distribution of 
transmission eigenvalues for low-energy chaotic scatter- 
ing in the four superconducting ensembles. From this 
distribution all moments of the thermal conductance of 
an Andreev quantum dot can be calculated. In the limit 
of a large number of scattering channels the phenomena 
of weak (anti)-localization and mesoscopic fluctuations 
are analogous to those for the electrical conductance in 
the non-superconducting ensembles. The opposite single- 
channel limit, however, shows striking differences. Most 
notably, in the absence of time-reversal symmetry, the 
thermal conductance distribution is either peaked or sup- 
pressed at minimal and maximal conductance, while the 
corresponding distribution of the electrical conductance 
is completely uniform. 

While Andreev quantum dots with multiple scattering 
channels can be realized in a two-dimensional electron 
gas with s-wave superconductors, the single-channel limit 
is out of reach in these systems in the absence of spin- 
rotation symmetry because of a fermion doubling prob- 
lem. We have shown how Majorana modes at the inter- 
face between different topological phases can be used to 
overcome this problem. 

In closing we point to the possibility to realize the four 
superconducting ensembles in graphene, where a strong 
proximity effect to s-wave superconductors has been 
demonstrated^ An Andreev quantum dot in graphene 
could be created using superconducting boundaries,— 
as in Fig. El Since spin-orbit coupling is ineffective in 
graphene, only the two ensembles which preserve spin- 
rotation symmetry (CQE and T-CQE) are accessible in 
principle. However, if intervalley scattering is sufficiently 
weak (on the time scale set by the dwell time in the quan- 
tum dot), then the sublattice degree of freedom can play 
the role of the electron spin. This pseudospin is strongly 
coupled to the orbit, so one can then access the two en- 
sembles with broken spin-rotation symmetry (CRE and 
T-CRE). 

It is an interesting question to ask whether the single- 
channel limit might be reachable in graphene. For the 
CQE and T-CQE we need strong intervalley scattering, 
to remove the valley degeneracy. For the T-CRE we need 
weak intervalley scattering, and could use the very same 
setup as in Fig. [5] One can then do without a topolog- 
ical insulator, because the helical Majorana modes exist 
also in graphene at the interface between superconduc- 
tors with a 7r phase differenced For the CRE, however, 
weak intervalley scattering is not enough. We would also 
need to convert the helical Majorana mode into a chiral 
mode, which we do not know how to achieve without a 
topological phase. 
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Appendix A: Calculation of the transmission 
eigenvalue distribution 

We briefly outline how to obtain the distribution ((5|) of 
the transmission eigenvalues from the invariant measure. 
(For a more detailed presentation of this type of calcu- 
lation we refer to a textbook^) One goes through the 
following steps. The polar decomposition of S provides 
us with a parametrization in terms of the transmission 
eigenvalues Ti and angular parameters pi . We express the 
invariant measure dfi(S) in terms of these parameters via 
the metric tensor m: d)j,(S) = VdetmJ^ dxi, where {x{\ 
denotes the full set of parameters {Ti,p{\ and m is de- 
fined by Tr (dS^dS) = J^ij mijdxidxj. Upon integration 
over the p^s we obtain the required distribution P({Tj}). 

Starting from the first step, the polar decomposition 
reads 



S = 



where the Ni x N2 matrix A has elements Ajk = y/TjSjk- 
Referring to Table ITT1 the transmission eigenvalues have a 
twofold electron-hole degeneracy in classes C and CI, as 
a direct consequence of the fact that the matrix elements 
can be represented by (real) quaternions. In addition, 
there is a twofold spin degeneracy because spin-rotation 
symmetry is preserved. In class Dili, the presence of 
time-reversal symmetry produces a twofold Kramers de- 
generacy of the transmission eigenvalues. (We focus on 
the situation where Ni and A 2 are even.) The unitary 
matrices U n and V n are orthogonal in classes D and Dili 
and symplectic in classes C and CI. They are independent 
in classes D and C. In class Dili one has = o^t/Jo^, 
while in class CI = U* . 

The following steps are straightforward, apart from one 
complication. In the polar decomposition, the set of Tj 's 
and the matrices U n and V n introduce more parameters 
than the number of independent degrees of freedom of 
the scattering matrix. The metric tensor, however, is 
defined through the derivatives of £ with respect to the 
set of its independent parameters. Keeping {T{\ in our 
parametrization, we define the angular parameters {pi} 
as independent combinations of the matrix elements of 
SU n = U^dUn and 5V n = V^dV n . In this way, the sub- 
sequent integration over these degrees of freedom does 
not involve dependencies on the TVs. The integration 
over these parameters thus only produces an irrelevant 
normalization constant and need not be carried out ex- 
plicitly. 




8 



1 S. Datta, Electronic Transport in Mesoscopic Systems, 
(Cambridge University Press, Cambridge, 1995). 

2 H. van Houten and C. W. J. Beenakker, Phys. Today 49 
(7), 22 (1996). 

3 Y. Imry and R. Landauer, Rev. Mod. Phys. 71, S306 
(1999). 

4 C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

5 T. Guar, A. Miiller-Groeling, and H. A. Weidenmiiller, 
Phys. Rep. 299, 189 (1998). 

6 Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 

7 F. J. Dyson, J. Math. Phys. 3, 140 (1962). 

8 M. L. Mehta, Random Matrices (Academic Press, New 
York, 1991). 

9 H. U. Baranger and P. A. Mello, Phys. Rev. Lett. 73, 142 
(1994). 

10 R. A. Jalabert, J.-L. Pichard, and C. W. J. Beenakker, 
Europhys. Lett. 27, 255 (1994). 

11 C. W. J. Beenakker. larXiv:0904.14"32l 

12 D. V. Savin and H.-J. Sommers, Phys. Rev. B 73, 081307 
(2006). 

13 B. A. Khoruzhenko, D. V. Savin, and H.-J. Sommers, Phys. 
Rev. B 80, 125301 (2009). 

14 V. A. Osipov and E. Kanzieper, Phys. Rev. Lett. 101, 
176804 (2008). 

15 A. Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 
(1997). 

16 P. Heinzner, A. H. Huckleberry, and M. R. Zirnbauer, 
Commun. Math. Phys. 257, 725 (2005). 

17 S. R yu, A. P. Schnyd er, A. Furusaki, and A. W. W. Lud- 
wig, larXiv:0912. 21571 

18 M. R. Zirnbauer, I arXiv: 100 1.07221 

19 M. Caselle, Mod. Phys. Lett. B 10, 681 (1996). 

20 The parameters /3 and 7 + 1 are the socalled root multiplic- 
ities of the symmetry class, see M. Caselle and U. Magnea, 
Phys. Rep. 394, 41 (2004). 

21 T. Senthil and M. P. A. Fisher, Phys. Rev. B 61, 9690 



(2000). 

22 A. Vishwanath, Phys. Rev. Lett. 87, 217004 (2001). 

23 J. T. Chalker, N. Read, V. Kagalovsky, B. Horovitz, Y. 
Avishai, and A. W. W. Ludwig, Phys. Rev. B 65, 012506 
(2002). 

24 P. W. Brouwer, A. Furusaki, I. A. Gruzberg, and C. Mudry, 
Phys. Rev. Lett. 85, 1064 (2000). 

25 M. Titov, P. W. Brouwer, A. Furusaki, and C. Mudry, 
Phys. Rev. B 63, 235318 (2001). 

26 G. E. Volovik, The Universe in a Helium Droplet (Claren- 
don Press, Oxford, 2003). 

27 X.-L. Qi and S.-C. Zhang, Ph ys. Today 63 (1) 33 (2010). 

28 M. Z. Hasan and C. L. Kane. larXiv:1002. 3895 

29 U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986). 

30 C. W. J. Beenakker, Lect. Notes Phys. 667, 131 (2005). 

31 The basis in which 5* is real is constructed by taking the 
combinations |e) + \h), i\e) — i\h) of the electron and hole 
states |e), \h) with the same spin and mode quantum num- 
bers. The corresponding creation and annihilation opera- 
tors are identical, so these basis states represent Majorana 
fermions. 

32 B. Beri, Phys. Rev. B 79, 214506 (2009). 

33 I. Serban, B. Beri, A. R. Akhmerov, and C. W. J. 
Beenakker, Phys. Rev. Lett. 104, 147001 (2010). 

34 C. Kallin and A. J. Berlinsky, J. Phys. Cond. Matter 21, 
164210 (2009). 

35 L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

36 H. B. Heersche, P. Jarillo-Herrero, J. B. Oostinga, L. M. K. 
Vandersypen, and A. F. Morpurgo, Nature 446, 56 (2007). 

37 J. Cserti, I. Hagymasi, and A. Kormanyos, Phys. Rev. B 
80, 073404 (2009). 

38 M. Titov, A. Ossipov, and C. W. J. Beenakker, Phys. Rev. 
B 75, 045417 (2007). 

39 P. J. Forrester, Log-Gases and Random Matrices (Prince- 
ton University Press, Princeton, 2010). 



